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Abstract 

In this paper we consider several instances of the k-center on a line problem where the 
goal is, given a set of points 5* in the plane and a parameter k> 1, to find k disks with centers 
on a line I such that their union covers S and the maximum radius of the disks is minimized. 
This problem is a constraint version of the well-known fc-center problem in which the centers 
are constrained to lie in a particular region such as a segment, a line, and a polygon. 

We first consider the simplest version of the problem where the line I is given in advance; 
we can solve this problem in 0{n\o^ n) time. We then investigate the cases where only the 
orientation of the line I is fixed and where the line I can be arbitrary. We can solve these 
problems in 0{n'^ log^ n) time and in 0(n* log^ n) expected time, respectively. For the last 
two problems, we present (1 + e)-approximation algorithms, which run in 0(i7ilog^ n) time 
and 0{-^n\o^ n) time, respectively. 

1 Introduction 

A common type of facility location or clustering problem is the fc-center problem, which is defined 
as follows: Given a set of n points in a metric space and a positive integer fc, find a set of 
k supply points such that the maximum distance between a point in S and its nearest supply 
point is minimized. For the cases of the L2 and Loo-metric the problem is usually referred to as 
the Euclidean and rectilinear fc-center problem respectively. The cluster centers (i.e., the supply 
points) can be seen as an approximation of the set P. Drezner [9] describes many variations of 
the facility location problem and their numerous applications. 

fc-center problems (as well as many other clustering problems) arise naturally in many applica- 
tions (like, e.g., in the shape simplification or in the construction of bounding volume hierarchies). 
They can be formulated as geometric optimization problems and, as such, they have been studied 
extensively in the field of computational geometry. Efficient polynomial-time algorithms have been 
found for the planar fc-center problem when fc is a small constant [6l [10] . Also, the rectilinear 
2-center problem can be solved in polynomial time, even when d is part of the input [2D|. However 
both the Euclidean and rectilinear (decision) problems are NP-hard, for d = 2 when fc is part of 
the input [HI [21] , while the Euclidean 2-center and rectilinear 3-center are NP-hard when d is 
part of the input [20]. Hwang et al. [16] gave an n'^'^^^-time algorithm for the Euclidean fc-center 
problem in the plane. Agarwal and Procopiuc [1] presented an n^^'' ' ^-time algorithm for solv- 
ing the fc-center problem in R'^ under L2 and Loo-nretrics, and a simple (1 -f- e)-approximation 
algorithm with running time 0(n log fc) -h (fc/e)'^('=' ). 
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Problems. In this paper, we consider a variation of the fc-ccnter problem with potential applica- 
tions in wireless network design: Given n sensors which are modeled as points, we want to locate 
k base stations (or servers), which arc modeled as centers of disks, for receiving the signal from 
the sensors. The sensors operate by self-contained batteries, but the servers should be connected 
to a power line, so they have to lie on a straight line which models the power line. Thus we define 
the k-center on a line problem as follows: Given a set S of n points (i.e., sensors) and an integer 
k > I, find k disks with centers (i.e., base stations) on a line such that the union of the disks 
covers S and the maximum radius of the disks is minimized. We will also study several variants 
of the problem depending on whether the orientation of the center line is fixed or not. 

Related results. The problem of finding the fc-centers on a line is a constrained version of the 
standard fc-center problems where the centers are constrained to lie in a specific region such as a 
line, a segment, and a polygonn. 

For fc = 1, Megiddo [ini showed that the 1-center constrained to be on a given line can be 
computed in 0{n) time; our problem is a direct extension of this for fc > 1. Hurtado et al. [15j 
considered the 1-center in a convex region instead of the line, and showed that the 1-center lying 
in the convex polygon of m vertices can be computed in 0{n + m) time. If the 1-center is restricted 
to lie in a set of simple polygons with a total of m edges, then Bose et al. [3] proved that it can 
be computed in 0{{n + m) logn) time, which is an improvement to the previous quadratic time 
algorithm [4]. They also showed how to preprocess the n points in 0{n\ogn) time and 0{n) space 
such that for any query segment the 1-center lying on the segment can be reported in O(logn) 
time. 

For fc > 2, there are a few results [8l[22l[25] which have been done mostly for the base station 
placement problem in wireless sensor network. Das et al. |22| studied more constrained problems 
where the k centers lie on a specific edge of a convex polygon and the corresponding disks cover 
all n vertices of the polygon. Their algorithm runs in 0(min(?T,^, nfclogn)) time. In [221 125j . they 
considered different constrained 2-center problems in which the centers lie in a convex polygon [25j 
or on a pair of specified edges of the convex polygon [22] , but the two disks should cover all points 
in the polygon, i.e., the convex polygon itself. 

Alt et al. [2] studied a similar problem where the goal is to minimize the sum of the radii 
(instead of minimizing the maximum of the radii); in fact, they minimized the cost function of 
the form r" for any fixed a > 1 under any fixed Lp metric, where represents the radius 
of each disk. In their model k is not part of the input. They presented an algorithm that runs 
in O(n'^logn) time to compute the optimal solution for a fixed line where c = 2 for a = 1 and 
3 < c < 4 for a > 1. For a horizontal line moving freely, they proved that the value of the 
optimal solution (as well as the location of the optimal horizontal line) cannot be expressed by 
radicals, and they gave an (1-1- £)-approximation algorithm requiring 0{^n^ log n) time. They also 
presented an (1 + e)-approximation algorithm for arbitrary line which takes 0{^n^ ^ogn) time. 

Our results. We investigate the fc-center problem for arbitrary fc > 1 where the k centers are 
collinear and should cover a set of n points. In Section [5] we first consider the simplest version 
of the problem where the line I is fixed (we will assume w.o.l.g. that the centers of the fc disks 
must lie on the x-axis). We can solve this problem in O(nlog^n) time. Next, we consider two 
problems: In Section [31 we look at the case where only the orientation of the line (. is fixed (we 
will assume w.l.o.g. that I is horizontal), and solve it in 0(n^ log^ n) time. In Section [4l we study 
the most general case where the line £ can have an arbitrary orientation. We solve this problem 
in 0{n^ log^ n) expected time. For these two problems, we also present (1 -I- e)-approximation 
algorithms, which run in 0{^n\o^ n) time and 0{-^n\o^ n) time, respectively. These results 
are the first on the fc-center problem on a line for fc > 2. Note here that fc is part of the input but 
not involved in the running time, and all algorithms also work for the problems under any fixed 
Lp metric. The results are summarized in Table [1] 
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fc-centcr problems on a line 


Exact algorithm 


(1 + ff)-approximation 


Centers on a fixed line 


0(nlog^ n)^ 




Centers on a line with fixed orientation 


0(n^ log'^ n) 


0((l/e)7ilog^ n) 


Centers on a line with arbitrary orientation 


0(71* log^ n) 


0((f/£-')?ilog'^7i) 



"O(nlogn) time for Li and Loo metric. 



Table f : Summary of the results. 

Preliminaries. Let S = {si, . . . , s„} be the set of input points, where Sj := {xi,yi). We assume 
that no two points in 5* have the same x-coordinate, no four points lie on the same circle and that 
they are sorted in non-decreasing x-order. Let £ be the line on which the disk centers have to lie. 
Denote by D{p, r) the disk centered at p with radius r; if the radius r is fixed, we denote this disk 
by D{p). 

2 fc-centers on a fixed line 

We first consider the simplest version of the problem where the line t is fixed. Without loss of 
generality, we assume that € is the x-axis. Let ymax max^ \yi\. Let r* be the minimum radius 
such that there exist k disks of that radius with centers on the x-axis and with union covering S. 
Then it is clear that r* > j/max- To find the minimum radius r* , we will perform a binary search, 
combining the results of [71 [TH] with an algorithm for the following decision problem: 

Given r > j/max, decide if there exist k disks of radius r with centers on the x-axis and with 
union covering S . 

2.1 The decision algorithm 

Let Xi{r) be the intersection of D{si,r) Di, i.e., Xi{r) = [ai{r) , bi{r)] is an interval on the x-axis 
with the two endpoints ai{r) < bi{r). A piercing set of the interval set X{r) = {Xi(ry\ is a set of 
points on £ such that every interval in X{r) contains at least one point in the piercing set. We call 
the minimum cardinality of a piercing set for X(r) the minimum piercing number k{r) of X{r). 
This number can be computed in 0(n log 7i) time by selecting a piercing point at the leftmost 
right endpoint, removing all the intervals containing the piercing point, and repeating the same 
process for the remaining intervals until all the intervals are removed, c.f. |171 [5]. 

It is clear that S can be covered by the union of k disks of radius r with centers on t if and 
only if the minimum piercing number k(r) is not larger than k. Thus we can answer the decision 
problem in 0(7ilog7i) time; if the right endpoints (and the left endpoints) of Xi{r) are sorted then 
we can even do it in 0(n) time. 

2.2 Computing the optimal radius 

If the piercing number at r = ymax is no more than k, we find that r* = 2/max- So, in what 
follows, we assume that r* > j/max, and that the intervals Xi[r) = [ai(r), 6i(r)] are sorted in non- 
decreasing order of their right endpoints bi{r). As r increases, ai{r) decreases and bi{r) increases. 
We can interpret ai{r') and bi{r') as the x-coordinates of intersection points of the hyperbola 
Xi : r^ — {x — XiY' = yf with the horizontal line r = r' in the (x, r)-plane, as illustrated in 
Figure [TJ Since no bisector between Si and Sj is horizontal, every bisector intersects the x-axis 
£ once and thus two hyperbolas Xi and Xj meet exactly once. If we consider Xi as two curve 
segments by cutting it at its lowest point, then we have 27i curves with the property that any 
two of them intersects at most once. We will denote the set of these curves (which are also called 
pseudoline segments) by X. 

As r increases, the minimum piercing number fc(r) for the interval system X{r) is non- 
increasing. Our goal is to find the radius r* such that k{r*) is the largest integer with k{r*) < k. 
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Figure 1: Intervals defined on a fixed line. 



Since k{r) can cliange (i.e., decrease) only when the relative order of the endpoints of the intervals 
in X{r) changes, we know that r* is a radius where two hyperbolas from X meet. A straight- 
forward way to find the radius r* is to compute all the radii at which two hyperbolas meet, sort 
them, and do a binary search over these radii with the decision algorithm of the previous subsec- 
tion; sorting all the 0{n^) candidate radii takes O(n^logn) time and the binary search requires 
0(nlog^ n) time because it invokes the decision algorithm O(logn) times. Thus the total running 
time is 0{ti^ logn). 

We can improve the running time to 0(ri, log^n) by using the result of [71 I18| which was 
developed for the slope-selection problem. In the slope-selection problem we are given a set of n 
points in the plane and an integer m, and we have to determine the line of the m-th smallest slope 
among all the lines spanned by two input points. In the dual setting, this problem is equivalent to 
finding the intersection point with the m-th smallest y-coordinate among all intersection points 
of the lines (dual to the input points). This problem can be solved in O(nlogn) deterministic 
time [71 118j. Since the curves in X we are dealing with are pseudoline segments, we can adopt 
these algorithms immediately. Thus we can do one step of the binary search in O(nlogn) time: 
choose the median radius r^^cd by the algorithm of [71 [18] , and decide if r* < r^cd by the decision 
algorithm of the previous section. As a result, we can find r* in 0(7ilog^ n) total time. The result 
is summarized in the following: 

Theorem 1 Let S be a set of n points in M? and /c € N. Let r* he the minimum radius such that 
there exist k disks of that radius with centers on the x-axis and with union covering S . Then r* 
and such disks can be computed in O(nlog^n) deterministic time. 

Remark. We can compute r* in O(nlogn) time if we use the Li or the Loo metric. In these 
cases the functions ai{r) and bi{r) are linear functions with a slope of —45° and -1-45°, respectively. 
There are no intersections among the functions ai{r) and no intersections among the functions 
bi{r). Consequently the sorted sequence of the ai(r)'s (resp., &i(r)'s) defined at r = ?/max remains 
the same for any r > ymax- So once we have the sorted sequences at the beginning of the algorithm, 
we do not need to sort the endpoints again when we make decision for a fixed radius during the 
binary search; this results in a linear-time decision algorithm for a fixed radius r. 

A radius r will be potentially tested during the binary search only if bi{r) ~ aj{r) for 
some i j- We can represent these radii implicitly in a doubly sorted matrix as follows: Let 
o-Tr{i)j 0-77(2)1 • ■ • I o,TT{n) be the increasing sequence of Oi(r) for some permutation tt at r = y,nax, and 
let &cr(i), ^(t(2)i • ■ • j ^cr(ri) be the decreasing sequence of &i(r) for some permutation tr at r = 2y,„ax- 
We define r{'!r{i),a{j)) as the radius r such that a7r(i)('') = ba{j}ir); if r(7r(i), ct(j)) < j/max then 
we set r(7r(i), (t(j)) := 2/max- It is easy to see that r{ir{i),a{j)) < r(Tr{i),a{l)) for any I > j, and 
r(7r(i), > r(7r(/i), a{j)) for any h < i. Consequently, if we put the r{Tr{i),a{j)) into an n x n 
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matrix with rows (ordered by tt) representing a^^j) and columns (ordered by a) representing 
we get a matrix in which each row and each column is totally ordered. With this doubly sorted 
matrix, we can search the m-th smallest entry in the matrix in 0{n) time fl3j (in particular we 
need to access only 0{n) matrix entries). Using the 0(n)-timc decision algorithm outlined above, 
we can perform the binary search for r* in O(nlogn) total time. 

Theorem 2 Let S be a set of n points in and fc G N. Let r* be the minimum radius such that 
there exist k Li-disks (resp. Lao-disks) of that radius, with centers on the x-axis and with union 
covering S. Then r* and such disks can be computed in O(nlogn) deterministic time. 

3 /c-centers on a line with fixed orientation 

We now consider the case where only the orientation of the line £ is fixed. Without loss of generality, 
we may assume that £ is horizontal. Let r* be the minimum radius such that there exist k disks 
of that radius with centers on a horizontal line and with union covering S. For i/q € M, denote the 
horizontal line y = j/o by £{yo)- 

Like in the previous case, we first develop an algorithm for the following decision problem, and 
then find the minimum radius r* by using this algorithm and the techniques in [71 [IS] : 

Given r > 0, decide if there exist k disks of radius r with centers on a horizontal line £{yr) and 
with union covering S. 

3.1 The decision algorithm 

Consider the disks D{si) :~ D{si,r) of radius r around the points Si g S. As the line £{y) moves 
from y = — oo to y = +oo, we maintain the minimum piercing number k{y) for the intervals 
Xi{y) := D{si) O £{y). Since all the disks must intersect £{y), it is actually sufficient to sweep t{y) 
from the topmost bottom point of the disks to the bottommost top point of the disks. Denote 
the y-coordinates of these two points by ys and yt, respectively. Then for any y from ys to yf, the 
interval Xi{y) of any point Si is non-empty. Set X{y) := {Xi{y) \ Si G S}. 

As the line £{y) moves from y = ys to y = yt, the minimum piercing number k(y) changes only 
when the relative order of the endpoints of two intervals changes. Thus the events correspond to 
intersections of circles bounding the D{si)'s, and the total number of the events is O(n^). If we 
have a dynamic data structure to maintain the minimum piercing number of intervals for each 
event in 0{T) (amortized) time, we can handle all the events in 0{n^T) time. 

Chan and Mahmood describe a data structure to maintain the minimum piercing number 
of a set of n intervals under insertions with 0(log7i) amortized time per insertion. Using this data 
structure, we can maintain k{y) during the sweep with O(logn) amortized time per event, so we 
can handle all the events in 0{n'^ logn) time. 

Maintaining the piercing number. For completeness, we briefly explain how the data struc- 
ture proposed by Chan and Mahmood [5] can be applied to our problem. For a fixed y, we consider 
the intervals Xi = [ai,bi\ for Si € S (omitting 'y' in the notation), together with two dummy in- 
tervals Xq = [— oo, — oo] and Xn+i = [oo,oo]. The greedy algorithm to compute the minimum 
piercing number (as described in the previous section) chooses the right endpoints of the intervals 
as piercing points. If bi and bj for i < j are two consecutive piercing points chosen by the greedy 
algorithm, then bj ~ min;{fe; | ai > bi}. Now we define next(Xi) Xj if bj ~ min/{6; | a; > bi}, 
and partition the intervals into groups each of which is a maximal set of intervals with a common 
"next" value. Chan and Mahmood observe that the intervals in a group appear in consecutive 
order when all {n + 2) intervals are sorted according to their right endpoints. Thus we can define 
a weighted tree T of intervals with the "next" relation such that a pair of consecutive vertices in 
a group is connected by an edge of weight 0, and if Xi is the last interval in its group, then the 
vertex for Xi is connected to the vertex for next{Xi) by an edge of weight 1. As a consequence, 
the minimum piercing number corresponds to the total weight of the path from Xq to Xn+i in T. 
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The tree T is implemented with a data structure for dynamic trees [26j that supports hnk, cut 
and path-length queries in O(logn) amortized time. The intervals in each group G are maintained 
in a balanced search tree Bq ordered by the right endpoints so that the operations such as search, 
split, and concatenation can be done in logarithmic time. All the intervals are also stored in a 
priority search tree Q ordered by the left endpoints with priorities defined by the right endpoints. 
This allows us to find the next interval for a given interval in O(logn) time. These priority search 
trees support insertion and deletion of intervals in O(logn) time each. For further details, refer 

to [g. 

As the sweep line £(y) moves from yg to j/t, we maintain the intervals on (.[y) in the dynamic 
tree T. First, at y = y^, we construct the data structures of all n intervals by simply inserting them 
in 0(n log n) time [5]. The line t{y) stops at each event Vs < y l£ yti which is the y-coordinate 
of the intersection of the circles bounding two disks D(si) and D(sj) for some i and j; we here 
assume that is in the left of Sj. When these two circles intersect, the relative order of the 
endpoints of Xi{y) — [ai{y),bi{y)] and Xj{y) = [aj{y),bj{y)] changes just after y. Then we have 
four different cases that the order change. Let y"*" := y + e and y~ := y — e for a small e > (see 
Figure ID). 




(a) (b) 

Figure 2: The four possible events, (a) Case (i) at the lower intersection point and Case (ii) at 
the upper intersection point, (b) Case (iii) at the left intersection point and Case (iv) at the right 
intersection point. 



Case (i) foi(y^) < cijiy^) and hi{y^) > aj{y^). Two intervals Xi and Xj start to overlap after 
y. If the next interval of Xi at y^ is not Xj, then we have nothing to do. Otherwise, Xj is no 
longer the next interval of Xi at y^ . By the definition of the group, Xi must be the last interval 
in the group G that Xi belongs to. First we find the new next interval Xi of Xi at y+ using Q, 
where / ^ j. Next we split the group G into two groups Gi := G \ {Xi} and G2 := {A'.J. Then, 
at y"*", the next interval for Gi is still Xj, and the next one for G2 is now Xi. So we merge G2 
into the group G' that Xi belongs to at y~ . For these changes, we update T, Bq, and Be- 

Case (ii) bi{y~) > aj{y~) and hi{y^) < aj{y~^). Two intervals Xi and Xj start to be disjoint 
after y. The interval Xj is not the next interval for Xi at y~ because they overlap, but Xj can 
be the next interval for Xi at y+. If it is indeed, then we know from the definition of the group 
that Xi must be the first interval of the group G of Xi at y~ . We split G into Gi := {Xi} and 
G2 := G \ {^i}- The next interval of Gi is now Xj, thus Gi is merged with the group G' whose 
next interval is Xj at y~. As did in Case (i), we update T, Bq, and Be to reflect the changes. 

Case (iii) ai{y~) > aj{y~) and ai(y+) < aj(y+). In this case the relative order of the left 
endpoints is reversed after y, i.e., ai{y~^) < aj{y~^). For this change, we have nothing to do. 

Case (iv) bi{y~) < bj{y~) and bi{y^) > bj{y~^). The relative order of the right endpoints is 
reversed after y, i.e., 6i(y+) > bj{y^). For the group G whose next interval is Xi at y~, the next 
interval must be changed from Xi to Xj. So we merge G into the group G' of Xj. For this change. 
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we update T, Bq, Bqi and Q properly. 

As a result, we can update all the data structures in 0(log?7) amortized time per event, so we 
can solve the decision problem for a fixed radius in 0{in? logn) time. 

3.2 Computing the optimal line and radius 

We now describe how to find the optimal line i{y*) and radius r* by a binary search that uses 
the decision algorithm for a fixed radius. To run a binary search, we need a discrete candidate set 
for the optimal radius r*, and we first study necessary conditions for the optimal line i{y*) and 
radius r* . 

The optimal line i{y*) must be immobilized, in the sense that if we translate it either upward 
or downward then the radius of the disks should be increased in order to cover S. So in every 
optimal configuration, there must be at least two points of S on the circles that immobilize £{y*), 
as illustrated in Figure [31 In the first two configurations i{y*) is fixed by two or three points 
on one bounding circle, and in the other configurations it is fixed by at most four points on two 
bounding circles. 




Figure 3: Five configurations of the optimal line and radius; symmetric ones are omitted. 

From these optimal configurations, we obtain 0{n^) candidates for the optimal radius as fol- 
lows: consider the discs Dij with s,; and Sj on its bounding circle. Their centers move along the 
bisector of Si and Sj, and thus in the (y,r)-plane we can define the function r = rij{y) of the 
radius of Dij with center of y-coordinate y. It is easy to see that is a unimodal function, i.e., 

is decreasing for y < {yi + yj)/'2 and increasing for y > {y.i + yj)/2. If we define rii{y) as the 
distance of Si to the line i{y), we get another unimodal function rii{y) (which is decreasing for 
y < yi and increasing for y > yi). 

We now split r^j into two monotone pieces: a decreasing piece r^j and an increasing piece r^, 
and set TZ' := {r'. | 1 < i < j < n}, 7^+ := {r+ | 1 < i < j < n}, and 7^ := 7^- U7^+. Then every 
optimal configuration illustrated in Figure [3] corresponds to an intersection point of two pieces, 
one from TiT and the other from TZ^; the first case (a) corresponds to the point r~- H for a pair 
of Si and Sj, the second case (b) corresponds to the point r,^ n for a triple of Si, sj, and s/j, 
and the remaining cases (c), (d), and (e) correspond to the points r^i H r^ , r^^ n r^, and r~j n r^j 
for points Si, Sj, Sh, and s;, respectively. Therefore, the optimal radius r* is the r-coordinate of 
one of the intersection points between TZ~ and TZ^. Both TZ~ and 7?.+ consist of 0{n^) curves, so 
the number of intersection points is O(n^). 

Consequently, we can determine r* by performing a binary search (discriminated by the decision 
algorithm) on the radii associated with the vertices of the arrangement of the curves in TZ. But 
the complexity of the arrangement is O(n^), so we find a way to compute the median of these radii 
without computing the arrangement explicitly. Since any two functions in TZ intersect at most 
once, we can pick the median of the radii from the arrangement in 0{n^ logn) time, by using the 
modified slope selection algorithm [71[TH1 as we did in Section With this median of the radii, 
we run the decision algorithm given in Section l3.II in 0{n^ logn) time. Thus we can perform a 
step in the binary search in 0{n^ logn) time, and r* can be found in 0{n^ log^ n) time. 

We just proved the following: 
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Theorem 3 Let S be a set of n points in and k E N. Let r* be the minimum radius such 
that there exist k disks of that radius with centers on a horizontal line and with union covering S. 
Then r* and such disks can be computed in 0{n^ log^ n) deterministic time. 

Remarks. This algorithm can be immediately applied to the problem in the ii-metric. For 
the Loo-metric, it is trivially solved in 0{n\ogn) time because the optimal horizontal line is the 
middle line between the lowest and highest points. 

4 /c-centers on a line with any orientation 

We now consider the case where the line £ can have arbitrary orientation and position. For k ~ 2, 
the problem is equivalent to the standard two-center problem [6lll0j. so we assume that k > 2. Let 
r* be the minimum radius such that there exist k disks of that radius with centers on a line and 
with union covering S. We will find r* in a similar way as before: we first design a deterministic 
algorithm for the following decision problem and then perform a randomized binary search over 
some candidate set of radii: 

Given r > 0, decide if there exist k disks of radius r with centers on a line and with union 
covering S. 

The decision algorithm for fixed r runs in 0(n^ logn) time, and the randomized algorithm for 
finding the optimal radius r* by a binary search takes 0(n^ log^ n) expected time. 

4.1 Decision algorithm 

Let i{6, h) denote the line with slope 5 and y-intercept /i, and let r > be fixed. We need to 
decide if there exists a line £((5, h) such that ^((5, h) intersects D{si,r) for all 1 < i < n and the 
minimum piercing number k{5, h) for the interval system {_D(s,;, r) C]t{5, ft,) | 1 < i < ri} is no more 
than k. Note that the set {D{si,r) : € 5*} is fixed and thus that if such a line ^(5, h) exists, then 
we can move, i.e., translate and rotate, the line without making any combinatorial change, until 
the line reaches one of the following configurations: (i) it contacts two disks tangentially, or (ii) 
it passes through an intersection of circles bounding two disks and is tangent to a third disk, or 
(iii) it passes through an intersection of two disks and another intersection of two different disks; 
refer to the dash- lined disks of Figure [3] (c) . (d) and (e) . respectively. Thus we need to check only 
the lines in such configurations for this decision problem. 

For the first case (i), there are only 0{n^) possible configurations-the bitangents of the disks 
D{si, r)-so we can check them all with the decision algorithm for the fixed center line. This takes 
O(n^logn) total time. For the second case (ii), we have O(n^) possible configurations, and again 
we can check them all; this takes 0{n'^ logn) total time. We now describe how to check the O(n^) 
possible configurations corresponding to the third case (iii). We fix an intersection point o of 
circles bounding two disks, and check all center lines passing through o. To this end we first sort 
all the intersections of the other circles in angular order around o; let pi , . . . , pm denote this sorted 
sequence (note that m = 0(n^)). We denote the line passing through o and pi by £{pi). In a 
next step we determine a maximal interval I < a < b < m such that for any a < i < b the line 
£{pi) intersects all n disks. This can be done in 0(n^ logn) time by a simple angular sweep. We 
now maintain the minimum piercing number for the intervals on the line i{pi) as we sweep it from 
i{pa) to i{pb)- We can maintain the corresponding data structures in O(logn) amortized time 
per event as described in Section [3l so we can check in O (71^ logn) time all center lines passing 
through a fixed intersection o. We run this algorithm for all intersection points of circles bounding 
two disks, thus we can check the O(n^) possible configurations corresponding to the third case in 
0(n^ logn) time. 
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4.2 Finding the optimal line and radius 

To find an optimal triple (6*, h*,r*), we use another form of binary searching which was applied 
to the slope selection problem by Shafer and Steiger [^3]. As before, we investigate all optimal 
configurations to get a discrete candidate set of optimal radius r*. Recall the optimal configura- 
tions for the center lines of fixed orientation, shown in Figure [3] In each configuration, we can 
slightly rotate the center line in clockwise direction without increasing the radius while the union 
keeps covering S. Thus we need more points on another disk to immobilize the line, and wc obtain 
optimal configurations defined by two or three disks of radius r*; in every optimal configuration 
with two disks, at least one of the two disks is such as (a) or (b) of Figure [3] that contains a 
diametral pair of points or a triple of points on its boundary. Every other configuration consists 
of three disks, each having one or two points on the boundary. 

From these optimal configurations, wc collect candidates for the optimal radius r* . We first 
consider the configurations defined by two disks. As noted above, every optimal configuration 
with two disks includes a diametral pair or a triple of points of S that lie on the boundary of an 
optimal disk. This pair or triple determines the radius of the optimal disk that it lies on. So we 
simply compute the radii from all pairs and triples of points and get 0{n^) candidates for the 
optimal radius. Every other configuration consists of three disks, each having one or two points 
on the boundary. We can interpret the radii in such configurations as intersections of triples of 
surfaces as follows: Let fij{6, h) be the radius of the disk whose center lies on the line £{5, h) of 
slope S and y-intercept h and whose boundary contains two points Si and Sj of S; if i is equal to 
j then fii{S,h) is the distance from Si to £{S,h). The graph of /y is a well-behaved low-degree 
algebraic surface in 3-dimcnsional {S, r)-space, and each triple of these surfaces provides only 
constant number of radius values. Let F be the set of these 0{n^) surfaces, and A{F) be their 
arrangement; the complexity of .4(F) is 0{n^) (see page 533 in ^^)- 

Now we may compute all 0{n^) radii from the vertices of the arrangement A{F) and perform 
a binary search on the union of two sets of the previously computed 0{n^) radii and of the just 
computed 0{nP) radii, which will take 0(n^ log n) time. However, we can do better if we adopt 
the randomization technique, as used in the slope selection problem by Shafer and Steiger |23j . 
Instead of computing all 0{nP) vertices and radii from the arrangement, we select uniformly at 
random n'^ triples of surfaces in F and compute the radii from the triples; each triple of the 
surfaces gives us constant number of vertices from A{F), so we get 0{n'^) radii in total. We sort 
these O(n^) radii together with the previously computed 0{n^) radii, in O(n^logn) time. Using 
the decision algorithm for a fixed radius (of the previous subsection), we now perform a binary 
search and determine two consecutive radii and r.^+i such that r* is between and r^+i. This 
takes 0(n''log^n) time in total. 

Since the vertices were picked randomly, the strip M^[r,;, r^+i] bounded by the two planes z := ri 
and z := r^+i contains only 0(71^) vertices of A{F) with high probability; in fact this is always 
guaranteed if we select V,{n'^\ogn) triples of the surfaces. So we can compute all the vertices 
in W[ri,ri+i] by a sweep-plane algorithm [24] in C>(n^ log n) time as follows: we first compute 
the intersection of the sweeping plane at z := ri with the surfaces in F. This intersection forms 
a two-dimensional arrangement of O(n^) quadratic closed curves and straight lines with O(n^) 
total complexity, so we can compute it in O(n^logri) time. We next construct the arrangement 
in VF[ri, ri_|_i] incrementally by sweeping from the intersection at z := towards z :— r^+i. As a 
result, we can compute the 0{n^) vertices (and the corresponding 0{n^) radn) in W[ri, r^+i]. The 
computation time depends on the complexities of the curve arrangements on two planes z := ri 
and z :— r^+i plus the complexity of the surface arrangement in the strip r^+i], which is 

0(n'*). Thus the time to identify all 0{n^) vertices lying in the strip is 0{n^ logn). 

As a final step, we perform again a binary search over these 0{n^) radii in W[ri, r^+i] to find 
r*, which takes 0{n^ log'^ n) time. We can find the optimal radius r* in 0(n^ log^ n) time with high 
probability, so this randomized algorithm to find the optimal radius takes 0(n^ log^ n) expected 
time. This result is summarized in the following: 

Theorem 4 Let S be a set of n points in and fc e N. Let r* be the minimum radius such that 
there exists a set of k disks of that radius with centers on a line and with union covering S. Then 
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r* and such disks can be computed in 0{n^ log^ n) expected time. 

5 Approximation algorithms 

We propose two approximation algorithms for the problem of computing fc-line centers for lines 
with fixed and arbitrary orientations, respectively. 

5.1 Fixed orientation 

We consider the case where the orientation of the line £ is given in advance. Without loss of 
generality, we may assume that i is horizontal. Fix an approximation parameter e > 0. Let h be 
the difference of the y-coordinates of the lowest point L and the highest point H in S. Clearly, the 
optimal radius r* is at least h/2. For 6 :— eh wc sample h/5 — 1/e lines of equal distance between 
L and H, solve the problem for the fixed lines in 0(n log^ n) time per line, and take the smallest 
radius r' among the solutions. Since the optimal line lies between two consecutive sampled lines, 
the radius r' is at most r* + 6, which is r* + eh < (1 + 2e)r*. This result is summarized in the 
following: 

Theorem 5 Let S be a set of n points in R^. Let k and e > 0. Let r* denote the minimum 
radius such that there exists a set of k disks with radius r* centered on some horizontal line that 
cover S. We can compute in time 0{-nlog^ n) a set of k disks with radius at most (1 + e)r* 
centered on some horizontal line that cover S . 

5.2 Arbitrary orientation 

In this section, we give approximation algorithms for the general case where the line i containing 
the fc-centers is arbitrary. We first give a constant-factor approximation algorithm, and we show 
how to use this result to get a fully polynomial-time approximation scheme. 

We denote by £* an optimal line, and we denote by r* the optimal radius of k disks centered 
at £* and containing S. 

Lemma 1 We can compute in time 0{n\ogn) a radius rc and k disks with radius rc and with 
collinear centers, such that these disks cover S and r* <rc < V2r*. 

Proof. The width of S is the minimum distance between two lines that contain S] we denote this 
width by w. We first compute a line such that the maximum distance between £j^ and any 
point of S is at most w/2; this computation can be done in time O(nlogn) by first computing 
the convex hull of S, and then by finding the width of this convex hull in linear time using the 
rotating calliper technique [27] . 
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Figure 4: Proof of Lemma [T] The black dots are the points in S, and the white dots are the points 
in Sw = piS). 

We denote by p the orthogonal projection to £^1- We project all the points in S and obtain 
a point set Sw = {p{s) \ s g S}. (See Figure HI) We solve the fc-center problem for Sw when 
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the centers are constrained to lie on f^, . It is the one-dimensional fc-center problem, which can be 
solved in 0{n) time after the points in Sw have been sorted [T2]. We denote by the optimal 
radius for this problem, and we denote by Cw sl set of k points such that S is contained in the 
union of the k disks with radius and center in Cw We have r^, < r*, because when we project 
the k disks of a solution of the original problem to iw, we obtain a set of segments with length 
2r* whose union contains Sw Wc now distinguish between two cases. 

First we assume that < w/2. Let s be a point in S. There exists c € Cw such that 
|cp(s)| < ryj. As |sp(s)| < uj/2, it implies that \cs\ < w/^/2. We have just proved that S is 
contained in the union of the disks centered at with radius w/^/2. As we noticed earlier that 
r* > w/2, we conclude that the set of disks centered at Cw with radius w/^/2 is a V2- factor 
approximation of the optimum. 

Now we prove the remaining case: we assume that > w/2. Let s be a point in S. There 
exists c S Cio such that cp{s) < r^. Since |sp(s)| < w/2, we get \cs\ < \/2r^. It follows that S is 
contained in the union of the disks centered at with radius \/2r^,, and we conclude using the 
fact that < H 

We now extend Lemma [T] into an approximation scheme. We first compute a radius Vc such 
that r* < Tc < \f2r* . The diameter d of 5* is the maximum distance between any two points in 
S. We compute a pair a,b S such that \ab\ = d; it can be done in 0(7ilogn) time in the same 
way as we computed the width. The lemma below handles the case where S is skinny. 

Lemma 2 We assume that 0<e<l. If d> 3rc, then we can compute in time 0{e~^nlog^ n) a 
set of k disks with collinear centers and with radius less than {l+£)r* , such that these disks cover 
S. 

Proof. Let A > be a constant, to be specified later. We scale S so that = 1, and we choose a 
coordinate frame such that a = (0, 0) and b = (1,0). When i E Z, we denote rm = iXerc- For each 
i such that —6rc < rm < 6rc, using the result of Sect ion [5T] wc compute a (l-|-e/3)-approximation 
for lines with slope rrii, and we return the covering with disks of smallest radius. As we consider 
only 0(l/e) slopes to^, it takes total time 0(e~^nlog^ n). We will now prove the correctness of 
this algorithm. 
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Figure 5: The black dots along £* represent C* , and the white dots represent C. 

Let £* be an optimal line; we write its equation y = m*x + y* . Consider the two axis-parallel 
squares centered at a and b with edge- length 2rc. (See Figure [5l) Since \ab\ = d = 1 and d > 3rc, 
a line with slope outside the interval [— 6rc, 6rc] cannot intersect both these squares, and thus this 
line is at distance more than rc from a or b. Since r* < rc, we know that i* intersects both of 
these squares and thus m* G [— 6rc,6rc]. So there exists j € Z such that mj € [— 6rc,6rc] and 
I m j — m* I < Xerc . 

Wc denote by C* C £* the centers of k disks with radius r* that cover 5*. All the points in 
S have a;-coordinates in [0, 1], so we can assume that the points in C* also have x-coordinatcs in 
[0, 1]. Let xi, . . . ,Xk € [0, 1] denote the x-coordinates of the points in C*. We have 

C* = {c* I c* = {x,,m*x^ + y*) and 1 < i < k.} 
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Wc introduce the point set C obtained by translating each point of C* vertically until it reaches 
the line with equation y ~ rrijX + y* . Hence we have 

C = {ci\ Ci ~ {xi,mjXi + y*) and 1 < * < k.} 

Notice that for all z, we have 

\c*c^\ < \m* ~mj\< Xerc < XeV2r*. 

So choosing A — l/(3-\/2), we get that \c*Ci\ < er*/3 and thus S is covered by the disks with 
radius (1 + e/3)r* centered at C C £j. So if we consider the (1 + £/3)-approximate radius rj that 
our algorithm returned for lines with slope rrij , we have 

rj < (l + £/3)(l + e/3)r* < {l+eY, 

which completes our proof. H 

In the following lemma, we handle the remaining case where S is fat. 




Figure 6: Proof of Lemma [3l The crosses along £* represent C*, and the white dots represent C. 

Lemma 3 We assume that 0<e<l. If d < Sr^, then we can compute in time 0(e~^rt log^ n) a 
set of k disks with collinear centers and with radius less than {l + e)r*, such that these disks cover 
S. 

Proof. Let A > be a constant, to be specified later. For each integer i > such that 6i = iXe < n, 
using the result of Section |5.1[ we compute a (1 + e/3)-approximation for lines making an angle 
9i with horizontal. Among all these (1 + e/3)-approximate solution, we return one with minimum 
radius. As there are 0(l/e) angles 9i, it takes total time 0(£~^n log^ 7i). We will now prove the 
correctness of this algorithm. 

We denote by C* a set of centers in an exact solution to our problem, and we denote by 
9* £ [0, tt) the angle that the corresponding optimal line £* makes with horizontal. Since r* < d, 
there exists a point o G C* that is at distance at most d from a. We choose i such that \9i — 9*\ 
is minimized, hence \9i — 9*\ < Xe. We consider the line £i through o making the angle 9i with 
horizontal. We denote by (p the rotation around o with angle 9i — 9*; hence ip{£*) = ii- 

We now prove that the disks with radius (l + e/3)r* centered at f{C*) cover S. So let s denote 
a point in S. There exists a center c G C* such that \cs\ < r*. We have \oc\ < \oa\ + \as\ + \sc\ 
and thus \oc\ < 3d. It follows that |c(^(c)| < 3d\9i — 9*\ and thus |c(^(c)| < 3dXe. Then we have 
|iy9(c)s| < r* + 3dAe, and since d < 3rc and rc < V^r* , we get |</'(c)s| < (1 + 9V2Ae)r*. Choosing 
A — l/(27\/2), we get that the disks with radius (1 + e/3)r* centered at ip{C*) cover S. 

Consider the approximate A:-center Ci that was computed for lines making an angle 9i with 
horizontal, and the corresponding radius. We have < (1 + e/3)(l + e/3)r* < (1 + e)r* , which 
proves that Ci is a (1 + £)-approximate k center for lines with arbitrary direction. H 

Combining lemmas [2] and [3l we obtain the following fully polynomial-time approximation 
scheme: 
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Theorem 6 Let S be a set of n points in M? . Let fc G N and < e < 1. Let r* denote the 
minimum radius such that there exists a set of k disks with radius r* and collinear centers that 
cover S. We can compute in time 0(£~^nlog^n) a set of k disks with radius at most (1 + e)r* 
and collinear centers that cover S . 
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